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Analytical solution of the equation of motion for a rigid domain wall in a magnetic 

material with perpendicular anisotropy 

M. C. HickejE 

Francis Bitter Magnet Laboratory, Massachusetts Institute of Technology, 

150 Albany Street, Cambridge, MA 02139 USA. 

This paper reports the solution of the equation of motion for a domain waU in a magnetic material 
which exhibits high magneto-crystalhne anisotropy. Starting from the Landau-Lifschitz-Gilbert 
equation for field-induced motion, we solve the equation to give an analytical expression, which 
specifies the domain wall position as a function of time. Taking parameters from a Co/Pt multilayer 
system, we find good quantitative agreement between calculated and experimentally determined wall 
velocities, and show that high field uniform wall motion occurs when wall rigidity is assumed. 

PACS numbers: 75.60.Ch, 75.60.Jk and 75.70.Ak 



The area of domain wall spintronics is currently enjoy- 
ing its heyday, both as a fruitful discipline for investi- 
gating how conduction electrons impart angular momen- 
tum onto lattice magnetization spins [l| and from the 
point of view of industrial application. Dynamical stud- 
ies in domain wall transport [2] have led to their use as 
memory bits [3, 4] while domain walls also play a central 
role in magnetic logic devices [5J . Controlling nano-pillar 
magnetization with electron current [a| has been widely 
demonstrated and forms the basis for magnetic random 
access memory. 

Many studies on domain wall motion necessitate a full nu- 
merical treatment of the Landau-Lifschitz-Gilbert (LLG) 
equation together with a description of the total mag- 
netostatic energy. While the starting descriptions of 
the magnetostatic energy are well understood, the fi- 
nal numerical simulation often lacks the transparency of 
a purely analytical treatment. Domain wall motion in 
Permalloy thin films is richly complicated by a variety of 
topological structures which can be nearly energetically 
degenerate. Complications of domain wall distortion un- 
der field include the Walker breakdown effect and more 
generally, oscillatory motion, contraction and expansion 
of walls which are commensurate with the emission of 
spin waves. These effects are instabilities and the treat- 
ment of the wall as a singular object breaks down as 
the wall dissipates energy to the lattice. While permal- 
loy is an attractive material from the point of view of 
low magnetization switching fields and low anisotropy, 
this type of non-linear behavior is best avoided for re- 
producible shuttling of domain walls down a patterned 
magnetic wire. In this paper, we focus on the description 
of domain wall motion in a perpendicularly magnetized 
material (such as a Co/Pt multilayer). We show that, 
having assumed a rigid wall profile and negligible wall 
distortion (negligible spatial dependence of wall tilt an- 
gle), an analytical solution of the equation of motion of 
the wall under field comes out, and there are well defined 



limits where the domain wall motion is robustly linear. 
The assumption of negligible wall distortion is justified 
in these materials because the easy axis of the system 
is always perpendicular to the direction of motion. We 
begin with the LLG equation 

^^7(MxH,;^)-^(MxM) (1) 

where 7 is the gyromagnetic ratio defined as 7 = g^^ (g 
is the electronic g factor, and /i^ Bohr magnetron) and a 
is the Gilbert damping. We write the effective magnetic 
field in the system as follows : 
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E(j is the energy density which contains the exchange, 
uniaxial and magnetostatic external field energies as de- 
scribed by equation [3] In spherical coordinates it is writ- 
ten as: 

Ed = A[{yef + sin^ 6l(V0)2] - K cos^ 9 - ^qM • H (3) 

where K is easy axis anisotropy constant, A is the ex- 
change constant and /xq is the magnetic permeability of 
free space while and (j) are the spherical polar angles of 
the magnetization. 
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The magnetization (M=(Ma;,Mj^,Mz)) can be written in 
terms of the spherical polar angles (in a cartesian vec- 
tor basis) as M = Afs(sin0cos(/), sin0sin0, cos0), where 
(j) = (j){x,t) and 9 = 9{x,t) are the azimuthal and polar 
angles, respectively. We can write the time derivative of 
the magnetization in the basis vectors of spherical po- 
lar coordinates (cm, e^, e^). This is a more convenient 
coordinate basis, because the magnetic state of the sys- 
tem can be described by two scalar fields, representing 
the spherical polar angles, in the above set of equations. 
Further, only two coupled equations in cj) and 9 are re- 
quired to describe the magnetostatics and dynamics (see 
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FIG. 1: (Color Online) The wall position is plotted as a func- 
tion of time for applied fields between and 500 Oe. The 
Gilbert damping constant is fixed at q = 0.016. The arrow 
marks the increasing field. 



for example, Thiaville et al. 7]). Equation [T] now reads 



^ I M^ainO^ 

(5) 

From this matrix equation, we have a system of two cou- 
pled partial differential equations, which are first order in 
time. We can eliminate from the system of equations, 
and we then arrive at the following more simplified equa- 
tion describing the time evolution of the magnetization 
angle Q ; 




1 



M,(l 



7 



1 bE 



fiQ sm f 



aj SE 
^0 S6 



(6) 



We calculate the effective magnetic field (Equation [21) by 
means of variational calculus, in the following way : 
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where repeated indices are summed over and we have 
a similar equation for the azimuthal angle, (j). We now 
evaluate these expressions using the definition of the to- 
tal magnetostatic energy from Equation [3] and substitute 
these evaluated expressions into Equation [51 
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We now write down the magnetization of the wall, 
as a magnetostatic solution, and assume that the wall 
is rigid and undergoes no distortion (i.e. V0 = 0, 
V^(/) = 0). The magnetization for a Bloch wall in 
a material with perpendicular easy axis anisotropy is 



taken to be M = M,(0, 1/ cosh(2^^^), tanh(i^^^^)) = 
Ms (sin e cos 4>, sin 6* sin cj), cos (/>) , where x is the central co- 
ordinate of the wall magnetization and Q(t) is the po- 
sition of the center of the wall. We use the follow- 
ing parameterization for the magnetization angle ; 9 as 
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and insert this definition into the 



equation of motion given by Equation[51 and arrive at the 
following first order equation : 

-Q= .j,,\ hM,ga7 + 2^tanh(^)[-g+^] 
Ms[i- + a^j /io ^ ^ 

This equation can now be integrated, and an implicit 
solution for the wall position versus time is found to be : 



(l + /3u)(u + /3u^) = /3e 
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where u = e~^^~'5)/^ and the constants A and C are 
defined below in terms of the parameters of the magnetic 
material and to is an arbitrary constant. We can solve 
this equation above (whose left hand side is cubic in u) 
to find the solution in the explicit form Q(t) — F(A,C,t). 
The result of this inversion is as follows : 
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where y is given by the following relation : 
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v = 



3((-7i - V7i^+4aj/9 )/2)^ 
-((-71 - V7F+4^)/2)i 

The quantities ai and 71 are given by (7/3 — 4)/3/3^ and 
(72/?2 - 10 + 27/3 - /3e^^^^(*+*"))/(27/33), respectively 
while 13 = [A — C)/{A + C) and we define the con- 
stants A and C, as follows ; A = —{\a'yHapp)/{l + a^), 
C = {\2a-i/Ms{l -f a^)^la){-K + A/A^) and we choose 
the boundary condition dQ/dt(t=0)=0. 
The results of this analytical model are plotted in Figure 
[1] and we see two distinct regimes - a non-linear region 
for t < 60 ns and a linear regime which takes over at 
timescales greater than 60 ns for all field values. The 
values used here for the calculation are taken from a 
Co/Pt multi-layer material system [8,] with perpendicular 
anisotropy, as follows : q:=0.016, 7=2.2x10^A~^ ms~^, 
/io=47rx 10~^ N A~^, exchange constant for Co ; A=3x 
10""J m-\ M^=1.5 MA uT^ , K(=Ke//)=0.3xlO^ 
Jm^'^ and A ^ ^J A/ K=l{] nm. Note that the perpendic- 
ular anisotropy constant here K is an effective anisotropy 
constant which takes into account the effect of the thin 
film demagnetization field. Using these materials param- 
eters, the dynamic wall velocity (v=dQ/dt)) versus times 
at various applied fields (from to 500 Oe) is shown in 
Figure [21 (a) and this gives steady state wall velocities 
in the region 0-0.5 ms~^. The field direction is chosen 
so that reverse saturation of the magnetization occurs as 
the wall moves in the positive x direction. The steady- 
state (t > 60 ns) wall velocity is plotted in Figure [H (b) 
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FIG. 2: (Color Online) (a) Plot of the instantaneous velocity 
attained by the domain wall under motion by applied field 
at a fixed Gilbert damping parameter oi ot = 0.016. The 
plotted wall velocities here are for applied fields Oe to 500 
Oe, the arrow indicates the increase in field magnitude. The 
flat region of constant wall velocity is preceded by a critical 
region, (b) Plot of wall velocity as a function of field H at 
differing Gilbert damping constants showing the onset of wall 
propagation which occurs when the critical field is reached 
from saturation. 



as a function of applied fields at differing Gilbert damp- 
ing parameters. These results show that the wall begins 
to move once a critical field is reached and that the wall 
velocity has an exponential dependence on field. Fur- 
ther, we plot the wall velocity in the steady-state regime 
at an applied field of 500 Oe against the Gilbert damp- 
ing parameter a, as shown in Figure [31 Here we find a 
linear relationship for small a which corresponds to the 
models developed by Slonczewski 9.] and others [10, HH , 
whereby one takes the precessional regime of steady-state 
wall translation (post Walker breakdown) and writes the 
wall velocity as : v = ^_7 _i ff ~ ^XaH^ and this lin- 
ear expansion is valid for small a. For a — 0.3 and at 
\H\ ~ 500 Oe, we have a wall velocity of ■^ 5 ms^^. This 
is in reasonable agreement with recently published results 
[l2l | on field driven walls in Pt/Co(0.5 nm)/Pt thin film 
systems. That work reported experimental wall velocities 
of ~ 8-10 ms~^ at 500 Oe with a Gilbert damping con- 



stant of about 0.3, having established anisotropy energy 
density, exchange stiffness and saturated magnetization 
all identical to that which we have used to parameterize 
our analytical model, the results of which are plotted in 
Figure [3] and its inset. 




FIG. 3: (Color Online) Plot of wall velocity at \H\ = 500 
Oe as a function Gilbert damping constant. The linear trend 
(dashed line) corresponds to the precessional regime for small 
a. The inset shows the field dependent velocity at a range of 
Gilbert damping parameters. This calculation used magnetic 
parameters from the Pt/Co(0.5nm)/Pt multilayer system of 
Metajcas et al. 1131. 



This correspondence arises in the linear regime, where 
the wall translates uniformly and the models neglect 
pinning due to defects. The linear regime occurs after 
Walker breakdown and in the limit of a perfect wire and 
corresponds to the precessional regime. 
In conclusion, we have calculated an analytical solution 
of the equation of motion for a undistorted domain wall in 
a perpendicularly magnetized material. This solution is 
constructed using first principles arguments from energy 
minimum considerations and the trajectories of the wall 
are completely specified by material parameters. Under 
the assumption of wall rigidity, we have linear wall trans- 
lation above a critical threshold where the wall position 
is exponentially dependent upon time. The values for 
wall velocities in the linear regime are in good agree- 
ment with previous experiments on field-driven walls in 
Pt/Co(0.5 nm)/Pt thin films, and the wall velocity is lin- 
early dependent upon Gilbert damping corresponding to 
precessional motion for small Gilbert damping constant. 
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